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ABSTRACT 

We present results from numerical modeling of emerging flux regions on the 
solar surface. The modeling was carried out by means of 3D radiative MHD 
simulations of the rise of buoyant magnetic flux tubes through the convection 
zone and into the photosphere. Due to the strong stratification of the convection 
zone, the rise results in a lateral expansion of the tube into a magnetic sheet, 
which acts as a reservoir for small-scale flux emergence events at the scale of 
granulation. The interaction of the convective downflows and the rising magnetic 
flux undulates it to form serpentine field lines emerging into the photosphere. 
Observational characteristics including the pattern of emerging flux regions, the 
cancellation of surface flux and associated high speed downflows, the convective 
collapse of photospheric flux tubes, the appearance of anomalous darkenings, 
the formation of bright points and the possible existence of transient kilogauss 
horizontal fields are discussed in the context of new observations from the Hinode 
Solar Optical Telescope. Implications for the local helioseismology of emerging 
flux regions are also discussed. 

Subject headings: Sun: atmospheric motions — Sun: activity — Sun: granulation 
— Sun: magnetic fields — Sun: interior — Physical data and processes: MHD 



1. Introduction 



The most prominent magnetic structures on the solar surface are bipolar active re- 
gions. These magnetic complexes are comprised of a hierarchy of magnetic structures of 
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different sizes and exist over a range of timescales. With unprecedented spatial resolu- 
tion, time cadence a nd a stable point-spread-function, the Solar O ptical Telescope (SO T, 



Tsuneta et al.l 120081 ) onboard Hinode (formerly known as Solar-B, iKosugi et al. 



20071 ) is 



providing new insights into the fl ux emergence process. Using Sto kes profiles taken with 
the SOT SpectroPolarimeter (SP, iLites. Elmore, fc Streanderl l200ll ). various studies have 
been car ried out to examine magnet ic flux emergence a t scale s of the surface granulation. 
Whereas ICenteno et al.l (120071) andlOrozco Suarez et al.l (|2008r) studie d flux emergence in a 
quiet Sun environment, llshikawa et al.l (120081 ) and lOtsuji et al.l (120071 ) examined flux emer- 
gence in a plage region and near sunspot penumbrae, respectively. Despite the different 
regimes, all three studies point to the importance of gr anular conve c tive fl ows on the prop- 
erties of emerging flux, as was previously reported by iDe Pontieul ( 120021 ) using filtergram 
observations with the Swedish Solar Telescope (SST). 

In tandem with the improved quality of observational data available, numerical simu- 
lations of magnetic flux emergence have, in recent years, become increasing sophisticated in 
terms o f the physics included in the models. For instance. ICheung. Schiissler. fc Moreno-Insertis 
( j2007al ) performed 3D radiative MHD simulations to model the rise of buoyant magnetic 
flux tubes through the near-surface layers of the convection zone and the overlying photo- 
sphere. The observational diagnostics of flux emergence in granular convection from these 
models compare fav orably with the aforementioned studies performed with SOT data (e.g. 
Centeno et aDl2007l ). 



More recently, several groups have begun to model self-consistently the magnetic connec- 
tion between the near-surface convection zone layers all to way to the corona. For instance 
Tortosa et al. (in prep.) extended the simulations of ICheung. Schiissler. fc Moreno-Insertis 
(j2007al ) to include the evolution of the emerging magnetic field in the o verlying chromosphere 
and transition region. iMartmez-Sykora. Hansteen. fc Carlssonl (120071 ) carried simulations of 
the emergence of twisted flux tubes from the convection zone into the corona, including 
effects such as thermal conduction along magnetic field lines and non-LTE radiative terms. 
These numerical studies yield important new insights for behavior of emer ging flux in the 
chro mosphere and above. A l thoug h the magneto-convection simulations of lAbbettl (120071 ) 
and llsobe. Proctor, fc Weissl (120081 ) were not initialized with buoyant magnetic flux tubes, 
both studies show that convective flows interactin g with an ambient field n aturally lead to 
the creation of small-scale magnetic loops (see also Stein fc Nordlund 2009 ) , which emerges 
through granular upflows and interacts with the pre-existing magnetic fields in the atmo- 
sphere. 



In ICheung et al.l (j2007al ) , we presented a study of magnetic flux emergence in granular 
convection by means of numerical radiative MHD simulations with the MURaM code (IVogler et al 
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20051 ). The simulations in that study began with a buoyant magnetic flux tube embedded 



in the near-surface layers of the convection zone. The simulations were restricted to flux 
tubes with initials longitudinal fluxes up to 10 19 Mx. In this paper, we extend our previous 
study to flux tubes with a fl ux on the order of 10 20 Mx, comparable to the flux content of 
mid-sized ephemeral regions (jHagenaarll200ll ). 

The paper is structured as follows. In section El we describe the setup of the numerical 
experiments. In section 13.11 we discuss the subsurface evolution of the rising flux tube. 
In section 13.21 we compare the observational characteristics of our modeled emerging flux 
regions with observational results with particuar emphasis on Hinode SOT results. Finally, 
in section @J we discuss the implications of the this work for our understanding of emerging 
flux regions. 



2. Simulation setup 



We first obtained a relaxed 3D non-magnetic model of the near-surface layers of the 
convection zone by means of numerical simulation with the MURaM code, which solves the 
radiative MHD equations on a Cartesian grid. The simulations take into account the ef- 
fect of radiative energy exchange and changes in partial ionization in the equation-of-state 
(EOS). The code has been previo usly used to inv estigate various aspects of solar surface 



uding solar 



magnetism inc 
plage regions (IVogler 



2004 



' aculae (Keller et al J|2004rl. magneto-con vection in the quiet Sun, 
Vogler et all 120051 : iShelyag et al. l2007h and mixed pola rity re- 



gions ( Khomenko et al.l 12005 ). umbral magneto-convection ( Schussler fc Vogle 31200^. solar 

pores ( Cameron et al.ll2007l )7the photospheric reversed granulation (|cheung. Schussler. fc Moreno-Insertis 

2007b ). the solar surface dynamo ( Vogler fc Schussler 2007 : Schussler fc Voglerl 20081) and, as 

a pred ecessor to the present study, magnetic flux emergence in granular convection (jcheung. Schussler. fc M 

2007al ). All these simulations have been confined to depths below the photospheric base of 

1.8 Mm or less. In order to acco mmodate for the deeper layers co nsidered here, we use EOS 

tables from the OPA L project ([Rogers. Swenson. fc Iglesiasl 1 19961 ) for a solar gas mixture 

with abundances from Anders fc Grevesse ( 19891 ). 



The simulation domain has horizontal dimensions 24 x 18 Mm 2 and a height of 5.76 
Mm. The horizontal and vertical grid spacings are 25 km and 16 km respectively. The 
level z = corresponds to the mean geometrical height where the continuum optical depth 
at 500 nm is unity (i.e. r 50 o = 1) and is located approximately 300 km (2 pressure scale 
heights) below the top boundary. The choice of this location for the top boundary was 
motivated by the need to keep the time step sufficiently large (At is as low as 0.01 s due to 
the high Alfven speeds in strong field regions) for the simulations to progress. Above the 
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top boundary, the magnetic field is matched to a potential field configuration at each time 
step. Periodic boundary conditions are imposed at the side boundaries and open boundary 
conditions allo wing for smooth in- /ou t-flow are imposed at the top and bottom boundaries. 



Using the model convection zone/photosphere as the ambient non-magnetic state, an 
initially horizontal, twisted magnetic flux tube was introduced into the convection zone at a 
depth of z — —3.9 Mm. The longitudinal and transverse components of the magnetic field 
have the form 



where r G [0, 2R ] is the radial distance from the tube axis and R the characteristic tube 
radius. For the two (runs A and B) simulations presented in this paper, R = 600 km 
and the field strength at the tube axis is B = 14 kG. The longitudinal flux of the tube is 
$o = / B[dS = 0.98ttRIBo = 1.55 x 10 20 Mx. The dimensionless twist parameter specifies 
the relative strength of the longitudinal and transverse field components. Simulation runs A 
and B have twist parameters of A = 0.1 and A = 0.25 respectively. 

Although the initial field strength (B = 14 kG) of the tube is several times larger than 
the strongest fields observed at the solar surface, it is important to note that the near-surface 
convection zone and overlying atmosphere are strongly stratified. In fact, the initial depth 
of the tube is 7.7 pressure scale heights below the base of the photosphere (r 50 o = 1)- At 
this depth, the plasma-/? (ratio of gas and magnetic pressures) at the tube axis is f3 « 30. 

To keep the divergence of the total (Maxwell + gas pressure + viscous) stress ten- 
sor continuous throughout the domain, the initial thermodynamic state of plasma enclosed 
within the tube was modified. We imposed a sinusoidal profile for specific entropy distribu- 
tion within the tube so that s = 11.912* at x = and x = 24 Mm (lowest entropy at side 
boundaries) and s = 12.11?* at x — 12 Mm, where 1?* is the universal gas constant. The 
prescription of such a distribution for the specific entropy translates into a density perturba- 
tion for the tube which is buoyant at x — 12 Mm and anti-buoyant at the side boundaries. 
Thus the differential buoyancy along the flux tube induces it to develop into an f^-loop 
configuration. 



For details, see ICheung et al. 




Bi(r) 
B e (r) 



B exp(-r 2 /R 2 ), 



(1) 
(2) 
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3. Simulation Results 

In section 13.11 we first discuss aspects of the subsurface evolution of the buoyant flux 
tubes in the simulations. Thereafter, in section [3721 we discuss the photospheric diagnostics 
of our simulated emerging flux regions and their relation to observations. 



3.1. Subsurface evolution 



3.1.1. Horizontal expansion, deformation and fragmentation 

Figure [1] shows a time sequence of cross sections of the tube at x — 12 Mm in Run A. 
Between t = and t = 33 min, this section of the tube had risen a distance comparable 
to its initial diameter. Two effects are of particular interest. First of all, the tube tends 
towards a fragmentation into two counter-rotating vortex rolls (identifiable as magnetic con- 
centrations regions of particularly high field strength) with opposite signs of out-of-plane 
component of the vorticity. This effect is well-known from earlier numerical studies of the 
buoyant rise of untwisted and slightly twisted flux tubes and is explained by considera- 
tion of th e aerodynamic pressure difference across the tube interior and the external flow 



around it ( Schusslerlll9 79: 



Fan et allied : Icheung et alJbood )" 



Moreno-Insertis fc Emonetll 19961 : lEmonet fc Moreno-Insertisll 19981 : 



In addition, the rise of the tube is accompanied by a strong horizontal expansion, which 
is expected for the adiabatic rise of a fluid parcel (magnetic or otherwise) over multiple 
pressure scale heights. By the time it has reached the photosphere, its aspect ratio is such 
that it appears more like a magnetic sheet than a 'typical' flux tube. This result confirms 
the scenario sketched out in Fig. 4 of ISpruit. Title, fc van Ballegooijenl (119871 ). wherein they 
describe how a flux tube rising over a few pressure scale heights must expand laterally to 
become a flux sheet. This has potentially important implications for helioseismic attempts 
to detect emerging flux regions before they erupt onto the solar surface, namely that a flux 
sheet model is perhaps a more appropriate model than a flux tube model for helioseismic 
inversions. 

A direct consequence of the strong expansion of the tube during its ascent toward the 
photosphere is the progressive weakening of the field with decreasing depth (in the absence of 
variations along the tube axis and under the assumption of flux freezing, B oc g). In contrast 
to this, mass conservation entails that the amplitude of flow velocities in the convection 
zone increases toward the surface (the Mach number A4 ~ 0.1 at the interface between 
the convection zone and photosphere). These two factors act together in such a way that, 
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Fig. 1. — Time sequence of vertical cross-sections (at x = 12 Mm) through the rising 
magnetic flux tube in Run A (twist parameter A = 0.1). The rise of the flux tube over 
multiple pressure scale heights is accompanied by a horizontal expansion which leads to a 
sheet-like structure below the photosphere. 
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Fig. 2. — Vertical velocity pattern (greyscale) and magnetic field distribution showing the 
passage of the buoyant magnetic tube through a layer in the convection zone at a depth of 
2.3 Mm. Blue and red contours respectively delineate positive and negative polarity regions 
with vertical field strength B z > 1 kG and the overlaid vectors indicate the direction and 
strength of the horizontal components of the field. For reference, the three vectors at the 
lower-right corners of each panel have horizontal strengths of 1, 2 and 3 kG. 
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with decreasing depth, the influence of the convective flows on the field evolution becomes 
increasingly important. 

Figure [2] shows the passage of the rising tube through a layer of the convection zone 
at a depth of 2.3 Mm. The vertical flow field (grey-scale) shows convective cellular pat- 
terns not unlike the surface granulation. At this depth, however, the upflow cells have a 
typical length-scale of ~ 4 — 5 Mm (as opposed to ~ 1 Mm at the surface) and the down- 
flow network already shows sign s of fragmentation into separated downflow channels (cf. 
Benson. Stein, fc Nordhmdlbood ). At t = 24 min, the top of the rising tube begins to cross 
this layer. As the magnetic tube continues to rise, the upflow cell created by the rising tube 
expands and sweeps aside previously existing downflows. A comparison of the buoyancy and 
drag (by the downflows) f orces experienced by a rising tube shows that for a magnetic tube 



to rise against downflows ( jParkerlll979l : iMoreno-InsertisI 1 1 9831 : iFan et al.ll2003l ; ICheung et al. 
2007al ) . it must have an internal magnetic field strength 



B > 



7T 



H \ 1/2 



(3) 



where Cd ~ 0(1) is the drag coefficient, 71 Chandrasekhar's first adiabatic exponent, H p 
the local pressure scale height, R the tube radius, and B eq the equipartition field strength 
(i.e. the magnetic field strength such that the magnetic energy density is equal to the kinetic 
energy density of a downflow with speed v). At a depth of 2.3 Mm, typically 90% of the 
downflowing plasma has a speed of 2 km s _1 or less. Taking v = 2 km s _1 and a mean 
density of (g) = 1.7 x 10 -5 g cm -3 at this depth, we have B eq & 3 kG. In Fig. [21 we 
have seen that the emerging tube is associated with an upflow cell several Mm in diameter. 
Taking R = AH p (a conservative estimate) and taking 2Cd1i/k ~ 1, criterion ([3]) tells us 
that the emerging magnetic field must have a field strength B > 1.5 kG. The flux tube in 
our simulation exhibits field strengths of 2 — 3 kG at this depth. In accordance with the 
above estimate, it is able to overcome the convective downdrafts and emerge. 



3.1.2. Effect of twist on the rate of flux emergence 

Figure [3] illustrates the rise of the magnetic flux by displaying plots of the longitudinal 
magnetic flux crossing the vertical plane x — 12 Mm as functions of time for the two 
simulation runs (i.e. J B x dydz). Curves labeled z > — 3 indicate the amount of longitudinal 
flux above a depth of 3 Mm (the flux tube initially resides just below this depth), which is 
almost the same for the two simulation runs. Since the initial buoyancy of the flux tubes is 
the same (to first order in twist parameter A), this is not a surprising result. Toward shallower 
layers, the results from the two simulation runs diverge from each other systematically such 
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that the rate of transport of longitudinal flux is always higher in the more twisted case (run 
B). In particular, the rate of longitudinal flux emerging from the convection zone into the 
photosphere (z > 0) in run B is almost double that of run A (~ 10 18 Mx/min as opposed to 
- 5 x 10 17 Mx/min). 



This finding is consistent with the results of iMurray fc Hoodl (120081 ). By means of MHD 
simulations of the buoyant rise of magnetic flux tubes through an idealized stratified layer 
mimicking the convection zone and overlying atmosphere, they found that magnetic tubes 
with higher degrees of twist (and therefore greater magnetic tension) have higher rates of 
emergence into the atmospher e. In their case, they can cast the results in the context of a 
magnetic buoyancy instability ( jAchesonlll979l : lArchontis et al.ll2004l ; Murray et al.ll2006 ). we 
cannot apply the same criterion to our simulations because it hinges on the assumption of 
adiabatic perturbations from an initially stationary equilibrium. These assumptions, in the 
presence of large entropy changes near the photosphere and ambient convective motion, are 
clearly violated during the evolution of the magnetic structure in our models. 



3.1.3. Coherent subsurface roots and serpentine field lines of an emerging flux region 

Figure [H shows three synthetic magnetograms [z —5 Mm, z = — 2 Mm and z = Mm) 
from the emerging flux region in Run A (A = 0.1) at t = 62 min. Also shown are magnetic 
field lines. Deep down at z = —5 Mm, the field lines group together as rather coherent 
bundles to give a tidy bipolar structure in the magnetogram. Following the field lines higher 
up near the surface, the horizontal expansion of the rising tube separates the field lines and 
the granular dynamics undulates them to form a set of serpentine field lines. This result is 



cons istent with the det ailed observations of emerging flux regions by IStrous fc Zwaanl (119991 ) 



and lPariat et al.l (120041 ). who report that magnetic field lines in emerging flux regions emerge 
already undulated. Our result naturally identifies convective dynamics as the cause of the 
serpentine structure of field lines. 

For a simple £~Moop emerging through the surface, the ratio of the total unsigned mag- 
netic flux at the surface (^unsigned = J Z=Q \B z \dxdy) and the longitudinal flux that has actually 
emerged ($i ong = J z>0 B x dydz, measured through a vertical plane normal to the tube axis) 
is ^unsigned/^iong = 2. In the case of flux emergence with serpentine field lines, this ratio 
is higher. Fig. [5] shows time plots of both the longitudinal flux crossing x = 12 Mm above 
the convection zone (red line) and the total unsigned flux crossing the base of the photo- 
sphere (black line) for Run A (A = 0.1). Throughout the period of emergence activity, 
^unsigned /^long ~ 10, which indicates that emerging field lines are, on average, undulated 
5 times by the convective flow. Although the flux emergence rate in Run B (A = 0.25) is 
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Fig. 3. — Longitudinal magnetic flux crossing the x — 12 Mm vertical plane. The dashed 
and solid lines indicate values from the runs A (A = 0.1) and B (A = 0.25) respectively. The 
curves labeled z > — 3 indicates the amount of longitudinal flux above a depth of 3 Mm etc. 
The amount of longitudinal flux in the photosphere is indicated by curves labeled z > 0. 
Fluxes are given in units of the initial longitudinal flux of the tube, $o — 1-55 x 10 20 Mx. 




Fig. 4. — The magnetic configuration of the emerging flux region from its subsurface roots to 
the photosphere (run A, A = 0.1) at t — 51 min. The synthetic magnetograms (grey-shading 
of B z on horizontal planes) at different depths show the transition from a coherent bipole in 
the subsurface roots to the complex mixed-polarity pattern at the surface. The serpentine 
appearance of the field lines near the surface highlight the effect of granular dynamics on 
the morphology of emerging flux. This figure is also available as an MPEG animation. 
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higher, the ratio of the two fluxes at any given time is also ~ 10. This common value of 
^unsigned/^iong between the two runs is attributed to the fact that the spatial extents of the 
two emerging flux regions at the surface are similar. 



3.1.4- Vertical transport of magnetic helicity 

In the limit of ideal MHD, the injection of magnetic helicity into an infinite half space 
(i.e. {{x, y,z)\z > Zp}) from t he bottom boundary (z = zq) is described by the Poynting 



theorem ( IBerger fc Field 1 19841 ): 



dH 



2 j[(A p • B)v z - (A p • v)B z ]dxdy, (4) 



where v is the plasma velocity and A p is a unique vector potential satisfying the boundary 
and gauge conditions 

z.(VxA p ) = B z , (5) 
z-A p = 0, (6) 
V-A p = 0. (7) 

where z is the unit vector in the z-direction. The first term in Eq. describes the bodily 
transport of twisted magnetic structures by vertical flows through the base plane (emergence 
term). The second term describes the transport of helicity due to the braiding of magnetic 
field lines by horizontal shearing motion at the boundary (shearing term). 

I n acco rdance with our periodic side boundary conditions, we use Fourier transforms (cf. 



Chadl200ll ) to calculate A p for a given B z distribution. The cumulative magnetic helicity 
injected above a depth of 3 Mm and above the base of the photosphere (z = 0) are plotted 
in the upper and lower panels of Fig. [6l respectively. The values of helicity are normalized 
to $o (square of the original longitudinal magnetic flux of the tube). 

We find that within the convection zone, at z = — 3 Mm, the bulk of the vertical helicity 
transport is due to the emergence term, whereas near the surface, at z = 0, shearing term 
dominates to carry most of the helicity transport. Again, the reason for this change can 
be put into the context of the horizontal expansion experienced by the tube. In deeper 
layers below the surface, the ratio of the tube radius to the pressure scale height is small, 
so that as parts of the tube traverses the z = — 3 Mm plane, it experiences relatively little 
expansion. As we found in section 13.1.11 the diminishing pressure scale heights near the top 
of the convection zone causes a dramatic expansion of the tube. In turn, the horizontal flows 
associated with this expansion enhance the shearing contribution to helicity transport. 
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Fig. 5. — Comparison between the longitudinal flux emerged into the photosphere 
($ion g = J z>0 B x dydz) and the corresponding unsigned vertical flux at the photospheric 
base ($ uns ig n ed = J \B z (z = 0)\dxdy). 
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Fig. 6. — Cumulative magnetic helicity injected above z = — 3 Mm (upper panel) and z = 
(lower panel) for simulation Run A (twist parameter A = 0.1). Helicity injection by the 
emergence term is indicated by dotted lines, the shearing term is indicated by dashed lines 
and their sum is indicated by solid lines. 
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In light of the work by Manchester et al. fl2004h and iMagaral (120061 ) , who reported that 
the expansion of rising twisted magnetic tubes into the solar atmosphere drives systematic 
shear flows, it is expected that the shearing t erm becomes increasingly important towards the 
surface. In fact. lMagara fc Longcopd (120031 ) finds from their emerging flux simulations that 
although the shearing term is initially smaller than the emergence term (for the photospheric 
base), the shearing term persists longer and becomes the dominant contributor to helicity 
injection into the solar atmosphere. 

Strictly speaking, Eq. (Ill) applies to an infinite half-space with no side boundaries. So 
the use of Fourier transforms for the computation of A p is only consistent with the scenario of 
identical emerging flux regions periodically spaced apart. To mimick the scenario of a single 
emerging flux region embedded in field- free surroundings, one can numerically evaluate the 
Biot-Savart integral for A p over the l ower boundary of an i nfinite half-space by setting B = 



outside the simulation domain (e.g. iBereer fc Field! 1 1984c iBergerl Il999l ; IWelsch fc Longcope 
20031 : IPariat. Demoulin. fc Bergerl 120051 : Ueong fc Chad 120071 ) . Application of the two meth- 
ods for this study show that the Fourier transform method typically yields helicity fluxes 
(and hence cum ulative injected helicit ies) which are ~ 10 — 20% higher than the direct inte- 
gration method. Ueong fc Chad (120071 ) also compared helicity fluxes from both methods and 
found a comparable excess from the Fourier transform method. 



3.2. Observational diagnostics: comparison between model and Hinode Solar 

Optical Telescope observations 

3.2.1. A magnetic inversion layer above granular up flows 

During the period of intense emergence activity, plasma both above and below the 
r 5 oo = 1 surface is permeated by magnetic field. This is shown by the cross-sectional plots 
in Fig. In the top panel, we find that below r 50 o = 1 (in granular upflows), the magnetic 
field strength generally increases with depth. This is consistent with the quasi-adiabatic 
expansion of frozen-in magnetic field in expanding upflows. Immediately above the T500 = 1 
surface on top of granules, however, is a layer where the magnetic field strength exceeds 
that of the plasma below the surface. This magnetic inversion layer is coincident with an 
inversion layer in the mass density (central panel of Fig. [7]) where material just above r 50 o = 1 
is top-heavy (see Fig. [8] for a vertical profiles of \B\ and g through a typical upflowing region). 
The relative increase of \B\ across the boundary layer (~ 50%) is larger than the relative 
increase of g (^ 30%). 

In non-magnetic solar convection simulations, such a density inversion layer above gran- 
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y [Mm] 

Fig. 7. — Vertical cross-sections (at x — 12 Mm) of magnetic field strength |£>|, mass density 
g and specific entropy s (top to bottom, respectively) at in the near-surface layers of the 
emerging flux region in Run A (twist parameter A = 0.1). In all three panels, the white solid 
line near z = indicates the r 50 o = 1 surface. In the bottom panel, the red arrows indicate 
the fluid velocity in the plane. Immediately above this surface is a magnetic inversion layer. 
Vertical profiles of \B\ and g along the black dashed line are shown in Fig. 




Fig. 8. — Vertical profiles of \B\ (solid line) and g (dashed line) taken along the vertical 
dashed line in Fig. [3 clearly indicate the presence of the magentic inversion layer. 



- 15- 



ules is w ell-known (see Fig. 18 of Stein fc Nordlundl ( 1998 ) and Fig. 8 of Cheung. Schiissler. fc Moreno-Insert 



( l2007bl ) for plots of the density history of fluid parcels emerging into the photosphere) . The 
high density layer results from the abrupt radiative cooling experienced by plasma emerging 
onto the photosphere (see the sharp drop in specific entropy in the lower panel). 

By combining the mass continuity equation and the Lagrangian form of the ideal MHD 
induction equation, 

DB 



S(V >v) + {B- V)iT, (8) 





Dt 

we can derive Walen's equation for B/g: 

?\ ( u \ 

(9) 

This shows that the compression of plasma also leads to an intensification of the predom- 
inantly horizontal magnetic field emerging onto the solar surface. In the absence of the 
stretching of field lines, the r.h.s. of this equation vanishes and we have B/g = constant. 
From Fig. [8] however, we see that the ratio \B\/g increases with z (especially pronounced 
above z = 0). This shows that the stretching of magnetic field lines by the convective flow 
1) reduces the weakening of the field due to expansion in the convection zone and 2) leads 
to the stronger relative change of \B\ across the inversion boundary layer in comparison to 
that of g. 



3.2.2. Complex patterns in the surface magnetic field distribution 

One of the most striking features of the surface magnetic flux distribution in our sim- 
ulations is the complexity of the mixed-polarity patterns (see Fig. [9]). Instead of a pair of 
coherent flux concentrations of opposite polarity, the simulated emerging flux region (EFR) 
consists of a large number of flux concentrations residing in the intergranular network of 
downflows. Although there is a net bipolarity of the flux pattern (there are more positive 
then negative polarity concentrations on the r.h.s. and vice versa), there is also a considerable 
amount of mixed-polarity in the interior of the emerging flux region. This is an observa- 
tional consequence of the undulation experienced by field lines rising against the convective 
downdrafts (see section 13.1.31) . The net bipolarity of the region becomes more apparent with 
sufficient spatial averaging (lower resolution). 

Although the mixed-polarity pattern is complex, the serpentine nature of the field lines 
gives it a certain order, namely that the positive polarity concentrations appear at the 
left edge of granules and stream leftwards, whereas negative polarity concentrations appear 
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Time — 50 min Time — 66 min 




Fig. 9. — Continuum intensity images (at 500 nm) and synthetic magnetograms of the 
simulated emerging flux region in Run B (twist parameter A = 0.25). The vertical component 
of the magnetic field (B z ) is sampled from the r 50 o = 0.1 surface along vertical lines-of-sight. 
The interaction of the flux tube with granular convective motions is clearly imprinted onto 
the surface flux pattern, which is far from a simple, tidy pair of opposite polarity flux 
concentrations. Magnetic bright points in the intergranular lanes can be discerned in the 
snapshot at t = 66 min. This figure is also available as an MPEG animation. 
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at the right edge and stream rightwards. The apparent horizontal motion of individual 
magnetic elements from the simulated EFRs have speed of about 1 — 2 km s _1 , with magnetic 
elements near the fringe migrating at a faster pace than those close to the ce nter of EFR. 



This i s consistent with previ o us reports of moving magnetic features in EFRs (jStrous et al. 



p 
i 



19961 : iStrous fc Zwaan 119991 : iBernasconi et al.l 120021 ) . In such a scenario where opposite 



polarity flux concentrations are counter-streaming, the encounter of opposite polarity fields 
is common. The effects of flux cancellation within an emerging flux region are discussed in 
the following section. 



3.2.3. Surface flux cancellation as a source of supersonic downflows 

Encounters between flux concentrations of opposite polarities within an EFR have a 
number of interesting observational consequences. Firstly, magnetic diffusion between op- 
posite polarity fields and the retraction of inverted U-loops leads to a decrease of the total 
unsigned flux at the photosphere. Secondly, such flux removal sites may be locations of 
particularly strong, in some cases even supersonic, downflows. Fig. [TU1 shows a small region 
within the simulated EFR in Run B (A = 0.25). Between the pair of opposite polarity 
flux concentrations in the center of the field of view, we find a vertical downflow with a 
speed exceeding 6 km s _1 (red solid). In the very center of this downflow, the speed reaches 
values of up to 10 km s _1 . While such a downflow is supersonic (Ai s ~ 1.3), it is still sub- 
Alfvenic (A4a < !)• This suggests that Maxwell stresses in the magnetic field are primarily 
responsible for the acceleration of the supersonic downflow. 

Figure [TT] shows a 3D rendering of the flux cancellation site depicted in Fig. [101 The 
semi-transparent grey-scale surface shows the footpoint of the opposite polarity magnetic 
concentrations (B z at z — 0) at the photospheric base. Magnetic field lines through the 
flux concentrations are color-coded by the local values of v z , with reddish colors indicating 
downflows with v z < —6 km s _1 . The turquoise colored sheets are region of enhanced j 2 /g 
(we divide by g to fold out the effects of density stratification). Clearly, the regions of 
particularly strong downflows are located at the places where the curvature of the field lines, 
and hence the magnetic tension, is strongest. A quantitative examination of the forces in the 
downflow region reveals that: (1) the downwards directed gravitational acceleration (gg) of 
the plasma is roughly balanced by the upwards directed gas pressure gradient force (— Vp gas ) 
and (2) the vertical component of the Lorentz force {c~ l j x B) is directed downwards and 
has an amplitude several times that of both the pressure gradient and gravitational forces. 
Thus the Lorentz force is clearly identified as the cause of supersonic downflows at flux 
cancellations sites in the photosphere. 
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Fig. 10. — Synthetic magnetogram (B z sampled at r 50 o = 1) of the interior of the simulated 
emerging flux region (Run B, A = 0.25) at t = 60 min. Yellow and red contours show 
downflow regions with vertical speeds of 2 km s _1 (yellow contours) and 6 km s _1 (red 
contours) respectively. This snapshot shows an example of a supersonic downflow at a flux 
cancellation site. A 3D rendering of this region is shown in Fig. [TIJ 




Fig. 11. — 3D rendering of the flux cancellation site depicted in Fig. [TQl The magnetic field 
lines are colored according to the local value of v Z) with red indicating a vertical velocity 
v z < —6 km s _1 . The turquoise colored sheets are regions of enhanced j 2 /g (i.e. sites of 
enhanced magnetic gradients). 
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Fig. 12. — Hinode SOT observation of the emerging flux region leading to the formation of 
AR 10978. The left panel shows an SpectroPolarimeter (SP) scan of the EFR (sampled in 
the blue wing of the Fe 6301.5 line). The top right panel shows a zoom view of the flux 
cancellation region. The lower right panel shows the Stokes V/I profiles along the dashed 
diagonal line in the overlying panel. The redshift of the profiles associated with a supersonic 
downflow can clearly be seen at the position of the magnetic polarity reversal. 
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Detections of supersonic downflows (v z < —6 k m s 1 ) in the sola r photosphere have 



previ ously been reported inside a delta sunspot pair (IPillet et al.lll994l ) and in the penum 



bra (Idel Toro Iniesta. Bellot Rubio. fc Colladosl l200lh of an individual sunspot. More re- 



cently, using Hinode SP data, IShimizu et al.l (120081 ) reported that such supersonic down- 



flows in fact appear frequently in a wide range of environments in and outside of sunspots 
and active regions. They sketched a number of scenarios which may cause such supersonic 
downflows, one of which involves magnetic reconnection. 

In Fig. [T2J, we show an example of a supersonic downflow which is most probably a 
result of reconnection in or above the solar photosphere. The left panel of this figure shows 
a Hinode SP scan of the EFR associated with AR 10978. The top right panel shows a zoom 
view of a small region where there is a clear reversal of the longitudinal magnetic polarity. 
In the panel immediately below, we show the Stokes V/I profiles sampled along a line indi- 
cated by the diagonal white dashed line. At the polarity inversion line, there exists strongly 
redshifted Stokes V signal. To examine this downflow in further detail, we refer to the plots 
of the individual Stokes profiles in Fig. [131 The spatial locations at which the Stokes profiles 
were sampled are indicated by the diamond (negative polarity), square (polarity inversion 
line) and triangle (positive polarity) symbols. One of the interesting features of the profiles 
at the polarity inversion (red curves) is the W-shape of the Stokes V signal, which is con- 
sistent with the fact that the pixel is at the p olarity inversion line and thus its profile s have 



contributions from opposite polarity fields ([Pillet. Lites. Skumanich. fc Degenhardtl 11994 
del Toro Iniesta. Bellot Rubio. fc Collados 200 ll ). The presence of significant linear polar- 



ization indicates horizontal components of B. Finally, the strong redshift of the profiles 
(especially in Stokes I and V) indicate plasma with downflows speeds approaching ^10 km 
s _1 . All these features are consistent with the scenario of flux cancellation in our simulation 
(see. Figs. HO] and E]). 

Before proceeding to other observational diagnostics of flux emergence, it is important 
to clarify that although the magnetic loops shown in Fig. [TT] are the result of reconnec- 
tion between field lines originally pertaining to distinct flux concentrations, the simulations 
we have performed here do not allow for a detailed study of the reconnection process it- 
self. Due to the location of the top boundary [z = 300 km) above the photospheric base, 
the connectivity of the photospheric field over length scales I > 300 km is determined by 
the potential field boundary condition. Consider two opposite polarity magnetic concen- 
trations in the photosphere whose field are originally not mutually connected. As the sep- 
aration of the flux concentrations decreases, their fluxes are likely to become connected 
by virtue of the potential field condition even before they are so close that reconnection 
within the domain takes place to connect the two flux concentrations. In simulations of 
resistive flux emergence that do include (albeit in an idealized fashion) layers above the 
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Fig. 13. — Stokes profiles sampled at the three positions shown in Fig. [T21 The red lines 
indicate profiles at the polarity inversion line (red square symbol) and the blue and green 
lines indicate profiles from opposite polarity regions on either side of the inversion line. 
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photosphere (jlsobe. Tripathi. fc Archontisl 120071 ) . the reconnection process is modeled and 
tends to occur in the chromospheric layers. Both scenarios, however, give a similar end 
result: namely a set of inverted U-loops which represent new magnetic connections between 
the flux concentrations. It is the subsequent relaxation of these inverted U-loops, and not 
necessarily how they were formed, which is of primary concern as to the source of the su- 
personic downflows. Furthermore, an examination of the terms in the induction equation 
(—v z dB z /dz vs. magnetic diffusion term) for the scenario in Fig. [TOl reveals that flux removal 
is predominantly a result of the retraction of inverted U-loops. 



3.2.4- Convective collapse and the formation of kilogauss field and associated bright points 



Since convective collapse was first suggested as a possible mechanism for the intensi 
fication of photospheri c flux tubes to kG strengths (IParkerl Il978l : IWebb fc Robertsl Il978l : 
Spruit fc Zweibellll979l ). a large body of observational and theoretical work has been carried 
out to test and to refine the original theory. A particularly important refinement has been the 
consideration of the i nhibitive effects of later al radiative energy exchange between a flux tube 
and its surroundings (jVenkatakrishnanlll986l ) , which essentially violates t he adiabaticity con- 
dition usually assumed in the convective collapse process. Calculations by lRajaguru fc Hasan 
(120001 ) of convective collapse in axisymmetric flux tubes show that radiative energy exchange 
does indeed inhibit convective collapse in flux tubes with $ < 10 18 Mx. This result is consis- 
tent with the observationally determined flux-field strength relation by lSolanki et al.l (119961 ). 



The inhibitive effect of lateral radiative energy exchange on field intensification has also 
been studied in multi-dimensional numerical models of solar surface magneto-convection. In 
such dynamic MHD models, the magnetic flux concentrations never really begin from equilib- 
rium configurations and the collapse mechanism is often referred to as conv ective intensifica- 



tion (ISchiisslerll 19901 : iGrossmann-Doerth et al.lll998l : iSteiner et al 1 ll998h . In lCheung. Schiissler. fc Moreno 



(j2007al ). we examined the flux-field strength relation for our flux emergence simulation and 
found - in agreement with iRajaguru fc Hasanl (120001 ) - that convective intensification is in- 
deed inhibited by lateral radiative energy exchange for flux concentrations with $ < 10 18 
Mx. 

Here we present such an example of a flux concentration undergoing convective intensi- 
fication from the simulations. Fig. [TH shows two snapshots of a magnetic flux concentration 
resulting from the expulsion of emerge flux to the intergranular lanes in Run A. At t = 54.9 
min, we find a flux concentration trapped in a downflow vertex. In the absence of magnetic 
fields, such downflow regions are usually relatively dark. However, as the continuum inten- 
sity map shows, the flux concentration at this instant already has a brightness comparable 
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to neighboring granules. At constant geometrical height z = (photospheric base), the field 
strength is predominantly sub-kG. A couple minutes later, at t = 57.2 min, convective in- 
tensification has proceeded to further strengthen this flux concentration. In this 'collapsed' 
state, even at z = 0, most of the tube has field strength in excess of 1 kG. A consequence of 
the partial evacuation of the flux concentration is the creation of a local Wilson depression, 
which reaches down to a depth of z = —300 km at its center. The presence of such Wilson 
depressions also means that magnetically sensitive photospheric lines probe deeper into the 
photosphere (as indicated by the r 50 o = 0.1 surfaces in Fig. [TH which is close to where 
the response of the Fe 6301.5 and Fe 6302.5 lines are strongest). So, in addition to a true 
strengthening of the surface field during convective intensification, the partial evacuation of 
flux concentrations gives an additional apparent intensification (compare the red and black 
solid curves in the lower panel of Fig. (Hj). 

At t = 57.2 min, the magnetic flux concentration appears as a bright point, with 
continuum intensity hoo/(hoo) ~ 1-6. From the Eddington-Barbier relation, we know that 
for vertical lines of sight, the outgoing intensity I v is given by the source function S v at 
optical depth unity. So, even though the flux concentration, at z — 0, is relatively cool 
compared to surrounding plasma at the same geometrical height, its temperature at r 50 o = 1 
is close to 7,000 K. In contrast, at the geometrical depth of the Wilson depression [z = 
—30 km), the plasma i n the adjacent granules can be much hotter (T « 8, 000 — 9, 000 
K). iBellot Rubio et aL (2000) have modeled the structure of plage fields by inversion of 
observed Stokes profiles assuming a thin flux tube model and reported a similar result. So 
although the magnetic flux concentration is radiating profusely, its internal temperature 
(at t 5 oo = 1) can b e maintained above the average value of the quiet Sun by side-wall 
heating (ISpruitl 1 1 976l ) . 

In our simulation, flux concentrations undergoing convective intensification have down- 
flow velocities reaching down to v z = —6 km s _1 . Unlike the flux cancellation case described 
in section 13.2.31 (where the Lorentz force was responsible for the acceleration to supersonic 
speeds), the downflows in this scenario of flux concentration is caused by an imbalance of 
the vertical pressure gradient (— Vp gas ) against gravitation (gg), sufficient to accelerate the 
material to 6 km s _1 over one or two minutes. After the formation of the bright point, the 
magnetic concentration evolves continuously in response to the buffeting of the surround- 
ing granules, with its relative brightness staying above //(/) > 1. This type o f evol ution 
resembles the behavior of bright magnetic elements observed by iBerger fe Titld (119961 ) . 



Recently, iNagata et al.l (120081 ) reported on observations by Hinode SP of a photospheric 
flux concentration undergoing convective collapse. Using a two component Milne-Eddington 
inversion, they identified the intensification of the field strength from 500 G to 2, 000 G. 
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Time = 54.9 min Time = 5T.2 min 
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Fig. 14. — Convective intensification as the cause of bright point formation. Top row: con- 
tinuum intensity (at 500 nm). Second row: surface synthetic magnetogram (sampled at the 
7"500 = 0.1 surface). Third row: vertical cuts of the temperature (grey-scale) and velocity field 
(vectors) along the yellow dashed lines, with the elevation of the r 50 o =0.1 surface indicated 
as red curves. Fourth row: plots of \B\ along the yellow dashed line including (red curves, 
sampled at r 50 o = 0.1) and excluding (black curves, sampled at z = km) the effects of the 
local Wilson depression. 
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The downflow associated with the collap se grew to 6 km s -1 over a time period of 150s. A 
similar scenario has also been inferred by iBellot Rubio et al.l (120011 ) from the time evolution 
of Stokes profiles of infrared lines. All these observed parameters of the collapse process are 
in good agreement with the example we presented above from our simulation. 



3.2.5. Anomalous transient darkenings 



The observational study of small-scale flux emergence events in an EFR by lStrous fc Zwaan 

10 



( 119991 ) indicate that transient darkenings (a few Mm in length and with life times of 
min) are associated with the emergence of small-scale flux bundles emerging into the solar 
atmosphere. Their statistical study showed that the transient darkenings are aligned with 
upflows on the order of 0.5 to 1 km s _1 and that the lanes were usually flanked by the 
formation of bright grains at the ends of the darkening. Based on such diagnostics, they 
conclude that the transient darkenings correspond to the crests of rising loops, whereas the 
flanking bright points corresponds to the photospheric footpoint of such loops. Furthermore, 
they suggest that the 'bright grains appear while the magnetic field in the footpoints is 
concentrating by convective collapse'. 



The observations by IStrous fc Zwaanl (119991 ) are, to a large extent, confirmed by the 
high-resolution observations from Hinode SOT. For instance, Fig. [T5l shows fractional circular 
polarization (NFI 6302 A Stokes V/I, top row) and G-band intensity (BFI at 4305 A) maps 
of an EFR near disk center (the coordinates are in arcsec from disk center) . The data 
was reduced using the FG_PREP routine in Solarsoft ( IFreeland fc Handyill9981 ) and image 
alignment was performed by manually comparing the G-band and Stokes I intensity images. 

In all three G-band images, one can find darkenings (2 — 3 arcsec long, marked by 
yellow ellipses) in the EFR which are thicker and darker than typical intergranular lanes 
(though the dark feature at (x, y) = (— 34", — 12") in the last frame is actually a pore). 
Upon close inspection, one also finds that bright grains typically appear at the ends of such 
transient darkenings. Analysis of the dynamics in those regions in the simulation reveals 
that convective intensification (see section 13.2.41) is indeed the m echanism responsi ble for the 



formation of bright grains, thereby confirming the hypothesis of IStrous fc Zwaan 



As reported bv Icheung. Schiissler. fc Moreno-Insertis ( 2007a ). the rise of small-scale 
flux tubes, given sufficient twist, are also able to create disturbances in the granulation 
pattern. For instance, a tube with 10 19 Mx flux and a twist parameter of A = 0.5 remains 
sufficiently coherent during its rise through granular convection that its horizontal expansion 
at the photospheric base modifies the local granulation pattern and leads to the transient 



-26 - 



2007-03-08T04:1 6:01 



Z007-03-08T04:37:02 



2007-03-08T04:56:02 



III 


, I , 


, , , 


I , i i. i I , , i i • I i i i i I , . i ! 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 ! I ! II I 1 1 1 ! 1 1 1 1 1 1 1 1 



08T04:56:33 




0.04 
0.02 

5 

I o.oo 8 

o 
co 

-0.02 

-0.04 
1 .6 

1 .5 

1 .3 A 

v 

1 .2 x) 

c 
o 

1 .0 O 

0.8 

0.7 



-50 -4£ 



Fig. 15. — Hinode SOT maps of fractional circular polarization (top, NFI Stokes V/I) and 
G-band intensity (bottom, BFI) of an ephemeral region during intense emergence activity. 
There exists several mixed-polarity patches within the center of the emerging flux region 
(EFR), an imprint of the granulation dynamics on the emergence process. Observational 
signatures of flux emergence in the G-band include elongated granules along the direction 
connecting the two main patches of opposite polarity flux, transient darkenings (marked by 
yellow ellipses) and bright point formation. 
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appearance of dark upflows. A general disturbance of the granulation pattern in our simu- 
lated EFR can also be seen in Fig. [91 In the intensity image on the right, we seen that within 
the EFR, the granular structures tend to be larger and more elongated than the 'typical' 
quiet Sun granule. 

Similar observational features in emerging flux regions can also be found in our simulated 
EFR. Fig. [T6l shows snapshots of the vertical component of the magnetic field {B Z) sampled 
at T500 = 0.1) and continuum intensity distribution near a darkening within the EFR. The 
darkening is predominantly associated with upflowing plasma (v z up to 1 km s _1 ) whereas the 
flanking bright grains at the ends of the darkening are associated with the footpoints of the 
rising flux bundle. In our EFR simulation, however, there also exist transient darkenings (of 
comparable spatial extent and lifetime) which are predominantly associated with downflows 
having speeds of about 0.5 — 1 km s _1 . In addition, even those darkenings which are initially 
associated with upflows eventually (within 5—10 minutes) become associated with downflows. 
This transition reflects that fact that although low entropy material can overshoot into the 
photosphere, it eventually overturns and forms new downflow lanes in the granulation pattern 
(mass conservation entails that most upflowing plasma must overturn within a few pressure 
scale heights). 



3.2.6. Transient kilogauss horizontal fields 



Measurements of horizontal photospheric fields outside of sunspot penumbrae have thus 
far resulted in sub-kG field strengths. For instance, ground-based observations of the quiet 



Sun with the Advanced Stokes Polarimeter 



wit h field strengths up to 600 G (ILites et al 



(ASP) revealed internetwork horizontal fields 



19961 ) . Similarly, ASP observations of EFRs 



by ILites. Skumanich. fc Martinez Pilletl (119981 ) and iKubo. Shimizu. fc Litesl (120031 ) led the 
authors to conclude that flux emerges initially as horizontal structures with sub-kG field 
strengths, and that only after the emerged fields become vertical do super-kG fields exist. 

Recent observation s of ubiquitous horizontal fields pervading the photospheric sur- 
face ( IHarvey et all 120071 : lOrozco Suarez et all 120071 : ILites et al.ll2008l ) have revealed average 
horizontal fields on the order of tens of G. The lack of a dependence of the mean horizontal 
field strength with solar latitude is strongly suggestive of the existence o f a surface dynamo, 
which has in fact been numerically modeled with the MURaM code ( IVogler fc Schussler 
2007 ). In any case, ubiquitous horizontal fields resulting from the surface dynamo are lim- 
ited to field strengths much weaker than 1 kG. In the following, we report that transient 
super-kG horizontal fields can exist within emerging flux regions. 
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Fig. 16. — Example of a transient darkening in the simulated emerging flux region (Run 
A, A = 0.1). The upper and lower panels respectively show the vertical field component 
(sampled at r 50 o = 0.1) and emergent continuum intensity at 500 nm. The transient darken- 
ing corresponds to the crest of a granular-scale Q loop emerging at the photosphere (yellow 
contours indicate magnetic upflows with //(/) < 0.7 and |£>hor| > 100 G an d v z — 0.5 km 
s _1 ). The footpoints of the loop correspond to the opposite polarity flux concentrations at 
ends of the darkening (the green and red contours, respectively, enclose positive and negative 
flux regions with \B Z \ > 700 G). Although these footpoints are located in downflows, they 
are relatively bright (1/(1) = 1.4 — 1.6). 
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Fig. 17. — Converging horizontal flows from neighboring granules can intensify surface hor- 
izontal fields to field strengths above 1 kG. The vertical magnetic field (B Z) greyscale) 
is overplotted with the horizontal flow field (indicated by yellow arrows) and contours of 
|£>hor| = 700 G and 1200 G. All quantities were sampled from the r 50 o =0.1 surface. 
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Figure [T71 shows the surface vertical field strength [B Z1 grey-scale), horizontal velocity 
field (arrows) and horizontal field strength (red contours indicate |£>hor| = ^00 G and 1200 
G) within a portion of the simulated EFR (run B) at t = 60 min. In this example, we find a 
couple of strong concentrations of predominantly horizontal field at the edge of neighboring 
granulation cells. The field strength in these concentrations can reach up to 1.5 kG, which is 
comparable to the stro ng vertical fields inferred from model inve r sions of Stokes profiles ob- 



serve d in the quiet Sun (jLinlll995l : ISanchez Almeida fc Litesll2000l : iDommguez Cerdena et al 



20031 ). Taking the mean density at the photospheric base to be (g) = 2.6 x 10 -7 g cm -3 , and 
the surface flow speed of granular motion to be v = 4 km s _1 , the photospheric equipartition 
field strength is B eq = (^(g)) 1 ^ = 700 G. This means that a field strength of 1.2 kG 
significantly higher than the equipartition value. 

The mechanism for the intensification of horizontal fields to super-kG strengths is related 
to the convective intensification process that occurs for vertical flux concentrations in the 
sense that mass drainage from the loop is important. Consider a granulation-scale fMoop at 
the surface, with the footpoints of the loop located in the intergranular lanes. As a result of 
mass drainage through the footpoints of the loop, the horizontal portions of the loop become 
evacuated (g/g ~ 0.5 — 0.8) and allow two effects to work together to intensify the field. First, 
the decrease in the internal gas pressure allows the colliding outflows from adjacent granules 
to compression the fluid. Second, as the colliding flows encounter each other and are deflected 
to the side, the horizontal portions of the small-scale loop are stretched. These two effects 
are basically the first and second terms of Eq. (|BD respectively. As is the case for vertical 
flux concentrations, the partial evacuation of the loop results in a local Wilson depression, 
so that the r 50 o = 0.1 surface is closer to z = than z = 160 km. This allows photospheric 
lines to probe into slightly deeper layers and provides an apparent intensification. 

Although we have presented an example of how kG horizontal fields may be formed in 
an EFR, we do not claim that they are common. Rather, we are merely pointing our the 
possibility that, under some conditions, such structures may indeed exist and perhaps be 
detected. 



4. Discussion 

The results in this paper have interesting implications for our understanding of the 
subsurface evolution of emerging flux as well as their observational diagnostics near the solar 
photosphere. By comparing our model results with the new observational data from the 
Hinode Solar Optical Telescope, we are able to discern the physical processes underlying 
many observed properties of emerging flux regions. 
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For instance, the mixed field morphology in the interior of EFRs is reproduced by the 
simulations, which r eveal the serpentine nature magnetic fiel d lines as they attempt to rise 
towards the surface (jStrous fc Zwaanlll999l : IPariat et al.ll2004l ). From this work, we find that 
the undulation is naturally explained by the interaction of convective downdrafts and the 
rising magnetic field lines. Although mixed-polarity field exists within an EFR, the pattern 
has a certain order to it, namely that flux of one polarity tends to stream towards one 



direc t ion, whereas the opposi te polarity flux migrates in the opposite direction (jstrous et al 



1999 : iBernasconi et al J 12002 ). At each of the two opposing edges of the the EFR, flux of 
one particular sign tends to coalesce. When sufficient flux comes together, the formation 
of solar pores and perhaps sunspots can occur. In a striking example of flux emergence 
observed by SOT (see Fig. [15]), solar pores formed within one hour after the beginning of 
flux emergence. Due to the diminishing time-step imposed by the CFL criterion (mainly from 
the high Alfven speed within strong field regions) , our simulations have not yet progressed to 
the stage of pore formation. However, we do find a trend in both simulations runs of surface 
flux coalescence, leading to some individual concentrations with up to 10 19 Mx. It is our aim 
in future studies to advance these (or similar) simulations to the stage tha t pores form as 



a nat ural consequence of flux emergence (as opposed to the simulations of ICameron et al 



20071 . which started with a coherent vertical flux tube as an initial condition). 



One interesting observational feature is the presence of supersonic downflows (v z < 
—6 km s _1 ) at some flux 'cancellation' sites, whic h can be found in bo th our models and 
the Hinode SP observations (see section 15.2.31 and IShimizu et all 120081 ). A possible cause 
of supersonic downflows at flux cancellation sites is revealed by the simulations to be a 
consequence of the relaxation of inverse-U loops (resulting from changes in connectivity). 
The magnetic tension in these loops is sufficiently strong to accelerate plasma to supersonic 
(M s > 1), but sub-Alfvenic (Ma < 1), speeds. For the example shown in Fig. [TUl the 
predominant cause of surface flux removal is actually the retraction of inverted U-loops (as 
opposed to magnetic diffusion). 



Following the observational study by iNagata et al.l (120081 ) in which they showed a 
clear example of convective collapse/intensification at work, we present further evidence 
of how this pro cess leads to the am plification of surface fields and to the formation of bright 
points (see also lVogler et al.ll2005l ). In section 13.2.41 we identified a clear case of this process 
happening to a magnetic flux concentration in our simulated EFR. 

In the context of local helioseismology, the results from this study suggest that in the 
near-surface layers, EFR regions are likely to be in the form of a magnetic flux sheet permeat- 
ing the layers in the convection zone immediately beneath what magnetograms reveal as the 
surface manifestations of EFRs. The spread of the tube into sheet-like structure is due to the 
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strong horizontal expansion it experiences in the top few pr essure-scale heights of the convec- 
tion z one. This behavior, which was originally predicted bv lSpruit. Title, fc van Ballegooijen 
( 119871 ) , can also be found in previous idealize d simulations of flux tubes emerging from the 
convection zone to the solar atmosphere (e.g. lMagarall200ll : lArchontis et al.ll2004l ). 



Another finding which may be relevant for the local helioseismology of EFRs is the 
presence of a magnetic inversion layer immediately above the photospheric base (i.e. above 
7500 = 1) of upflow granules (see section 13.2. II) . This layer extends throughout the horizontal 
extent of the simulated EFR and is present throughout the emergence process. It is associated 
with low entropy plasma which has recently emerged at the surface and experienced rapid 
radiative cooling. The result of the radiative cooling is an enhancement of the plasma density 
and magnetic field strength relative to the plasma immediately below. The thickness of the 
transition layer is comparable to the thermal boundary layer (a few tens of km) and the 
relative enhancements of g and \B\ are of order 10 -1 . Thus, the magnetic inversion layers 
above granules may practically act as sharp boundaries for propagating magneto-acoustic 
waves, the effects of which may need to be carefully considered. 
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